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Abstract 

The director configuration of dischnation fines in nematic fiquid crystals in 
the presence of an external magnetic field is evaluated. Our method is a com- 
bination of a polynomial expansion for the director and of further analytical 
approximations which are tested against a numerical shooting method. The 
results are particularly simple when the elastic constants are equal, but we 
discuss the general case of elastic anisotropy. The director field is continu- 
ous everywhere apart from a straight line segment whose length depends on 
the value of the magnetic field. This indicates the possibility of an elongated 
defect core for dischnation lines in nematics due to an external magnetic field. 
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I. INTRODUCTION 



Nematic liquid crystals are systems which are positionally disordered, but reveal a long- 
range orientational order This property is described on a mesoscopic level by a unit vector 
field n(r), which is called director. Due to the absence of permanent dipolar moments in 
nematics the director just indicates the orientation, but its has neither head nor tail. This 
particular feature yields very interesting defect structures in nematic liquid crystals. For 
instance, the director field shows line defects in three dimensions (or, equivalently, point 
defects in two dimensions), called disdinations, which have been studied and classified by 
topological methods Unlike in spin systems, disdinations of topological charge ±| are 

possible and stable in nematics. When an external magnetic field is applied perpendicular to 
such a disclination line, the resulting director configuration becomes even more interesting - 
it can be regarded as a domain wall filling a half-plane which terminates in the disclination 
line. Such walls with edges, known as planar solitons have been discussed for superfiuid 
Helium-3 by Mineyev and Volovik |^ . Whereas the qualitative behaviour of these soliton-like 
objects is well-established, a quantitative understanding of their structure can be obtained 
only from a thorough analysis of the underlying field theory. Its is the aim of our paper to 
perform such calculations. Our approach is based on a polynomial expansion of the director 
field. This method has been used previously both in relativistic field theories [Bipi and for 



the evaluation of domain wall dynamics in nematics [|T0[- It yields approximate analytical 
solutions for the director orientation. 

The paper is organized as follows. In Section II the director field equation for the 
planar solitons is derived. Section III develops a method for obtaining an approximate 
solution for the tilt angle of the director. Our technique is a combination of the polynomial 
expansion |]8|-|lO|] with further approximations which are tested by means of a numerical 
shooting method ||TT|. The discussion is performed within the framework of the Oseen- 
Zocher- Frank elasticity [Il2|-|T^. In Section IV we estimate the energy of the defect core of 
the planar solitons, and we minimize the total energy of the solitons in order to determine 
the length of the core. Section V contains concluding remarks. 



II. DIRECTOR FIELD EQUATION AND BOUNDARY CONDITIONS FOR 

PLANAR SOLITONS 

A. Free energy and field equation 

The geometry for planar solitons in nematic liquid crystals is drawn schematically in 
Figs. 1 (positive soliton) and 2 (negative soliton). The director field is essentially planar, 
perpendicular to a disclination line of strength ±| along the z direction of a Cartesian 
coordinate frame. Because the structure is independent on z, we restrict ourselves to the 
x-y plane {z = 0). Now we impose a magnetic field in the plane of the director along the 
X axis. Due to the magnetic anisotropy of the nematic the director tends to align along 
the magnetic field. However, the topological charge of the disclination has to be conserved. 
The resulting structure is a planar domain wall of Neel type 1 15,1^, 10 1 which ends in the 



disclination line Locally, close to the disclination, the director field preserves the defect 
structure. However, in a plane at a finite distance from the disclination line, which is given 
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by the half width of the planar Neel wall (Figs. 1 and 2), the director field is aligned 
parallel to the external magnetic field [0,^1. 

Due to the translational symmetry along the z axis it is sufficient to perform the calcu- 
lations in two dimensions only. The director orientation is then completely determined by 
the tilt angle field ^{x,y), which is measured with respect to the direction of the magnetic 
field H [x axis), 

n = cos y) a; + sin y) H = HqX. (1) 

We look for static director configurations, hence $ does not depend on time. 

The static director orientation inside the soliton corresponds to a configuration mini- 
mizing the total free energy F (per unit length in the z direction) which contains both the 
energy of the nematic phase F^em and the core energy of the disclination Fcore- (Within the 
defect core local phase transitions may occur.) The nematic energy -Fnem is the area integral 
of a free energy density J-'nem- This free energy density, in turn, consists of elastic contribu- 
tions (due to distortions of the director field) and of a magnetic part (taking into account 
the interaction of the nematic with the external magnetic field), hence ^nem = -^eiast + ^mag- 
The elastic free energy density follows from the Oseen-Zocher- Frank expression |T2|-Q. 

J^ciast = ^ Ku (div nf + ^ i^sg (n x curl nf . (2) 

In (^ Ku and K^^^ denote the elastic constants for splay and bend deformations in the 
nematic. Due to the restriction to planar director fields according to there are no twist 
deformations and the elastic constant K22 does not enter the calculations. 

The magnetic free energy density couples the director n to the magnetic field H via the 
anisotropy of the magnetic susceptibility (yUg means the magnetic field constant). 

^mag = -^/io Ax(n ■ i3")l (3) 

When inserting the ansatz for the planar director field (|1|) into (^ and (|^), we obtain the 
free energy density T^era of the nematic phase. 



J^nem = - (i^ll + K^z) ($' +$') + - (^^33 " ^ll) ($' " $') COS 2$ 

4 ^ A " 

+-(X33-i^ii)$x$, sin2<l>--/ioAx//o'(l + cos2$). (4) 

2 4 

In (^ and $y denote partial derivatives of the tilt angle with respect to the spatial 
coordinates. The energy of the defect core F^ore will be discussed separately in Section IV. 

The director configuration for the planar soliton, which minimizes the energy of the 
nematic phase, follows as a solution of the corresponding Euler-Lagrange equation 

nem 9^ aera q ( dj~ acni \ ^ 

Sui ~ dui ^ \d{d.jni) J 
The resulting equation for the tilt angle field $(a;, y) can be written in the following form 
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where 



X COS 2$) dy {^y COS 2$)] 
+K [dx {^y sin 2$) + dy (^^ sin 2$)] +K{^l-^l- 2$^$^) sin 2$ 
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sin 2$ =0, (6) 



i^33 - ^11 



Ku + K33 
is the elastic anisotropy. 



B. Boundary conditions 

The boundary conditions are an essential feature of the planar solitons. As discussed 
above, the defect structure is surrounded by a homogeneous director field and by a planar 
Neel wall. According to the choice of our Cartesian coordinate frame (Figs. 1 and 2) the tilt 
angle should be zero aX y = ±?/o, where Uq is the half width of the Neel wall. Additionally, 
it should glue smoothly to the homogeneous orientation. Thus the boundary conditions in 
the y direction (perpendicular to the magnetic field) are given by 

$(x, y = yo) = 0, $(x, y = -yo) = ±n, y = ±yo) = 0, (7) 

where ±7r is for the positive and negative soliton, respectively. 

In the X direction (parallel to the magnetic field), the director field at x < coincides 
with the planar Neel wall. For increasing x coordinate the domain wall structure is destroyed 
and the director field changes towards the homogeneous orientation, parallel to the magnetic 
field, which is reached at Xq. Hence, 

$(x = 0,y) = $Ncei(y), = xo, y) = 0. (8) 

It is important to note that the value of xq is yet unknown at this stage. 

The function $ncci(2/) describes the director inversion due to the planar Neel wall. It 
can be determined by solving the field equation (^ in one dimension. For x < there is no 
dependence on the coordinate x and the field equation is simplified to 



/io Ax 



1 - cos 2<^)%y + K<^1 sin 2$ - ^ " sin 2$ = 0. (9) 



Kn + K: 
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The center line of the cross section of the domain wall with the x-y plane coincides with 
the line x < 0, y = 0, with tilt angle $ = ±^ on it. Now, following the approach developed 
in our recent publication , we apply a polynomial expansion of the tilt angle up to third 



order in the distance y from the center line, 

2 4 2/0 4 \yoJ 

The different signs are valid for positive and negative solitons, respectively. Due to the choice 
of the coefficients in the expansion (p!oD , the boundary conditions (|^) are fulfilled. With the 
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approximate expression (0) for $Neei(|/) we can satisfy (|) up to terms proportional to the 
first power of y. This fixes the half width uq of the planar Neel wall, 



It is of the order of the magnetic coherence length. ( ]TT| ) and (|T^) are used in (^), which now 
provides the boundary conditions in the x direction. 

The approximate solution (^) could be improved by taking a higher order polynomial. 
If it is of the order y", then (j^) can be satisfied up to terms proportional to ?/"~^. In the 
present paper we shall restrict ourselves to cubic polynomials in which are sufficient to 
reveal our method of obtaining the approximate director field for the planar soliton. 



III. TILT ANGLE FIELD FOR PLANAR SOLITONS 

Our strategy for solving the non-linear partial differential equation (P) for the tilt angle 
y) proceeds in two steps. First we apply the polynomial expansion of the tilt angle field 
in the y coordinate. After separating the y dependence, we are left with a set of ordinary 
differential equations which is solved both numerically and, approximately, analytically. Of 
course, the polynomial expansion in y must satisfy the boundary conditions (|^. Therefore, 
up to third order (in congruence with the expansion for the Neel wall (PH])) it reads 

y) = (±^^ + C{x) y^ (y T y^f mod vr, for y > 0, y < resp. (12) 

The polynomial expansion (p!2D contains two unknown functions ^q{x) and C{x) that 
depend on the x coordinate. We can derive boundary conditions for them by inserting (|12|) 
into (P). This yields (for the positive solitons) 

^o{x = 0) = ^, %{x = xo) = 0, (13) 



2' 

T 

4y!' 



C{x = 0) = C{x = xo) = 0. (14) 



Our ansatz (0) is continuous everywhere apart from the x axis {y = 0). When crossing 
the X axis between x = and x = xq, a jump in the director orientation from +^q{x) to 
—$0(3;) occurs. This is connected to the physical singularity of the disclination line in the 
center of the defect. Most significantly, due to the infiuence of the external magnetic field 
the cross-section of the defect core is no more a point-like object in the x-y plane, but it 
is extended to a segment of a straight line of length xq. However, although the core of the 
defect is now strip-like (if we take into account the z direction), one can define its center 
line. It is located at x = Xd, where ^o{x = Xd) = f , which gives the largest jump (equal to 
^) in the director orientation at ?/ = 0. At x < there is no physical singularity, because 
$0 = — f is equivalent to $0 = +f • 

The discontinuity of (|T2p at y = refiects the fact that the continuum approach is no 



more valid close to the defect core, where strong gradients of the orientational order are 
apparent. On a molecular length scale around the core the mesoscopic director looses its 
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physical significance as the average molecular orientation. Remarkably, although when using 
the director approach we cannot determine the orientational order within the defect core, 
our investigation gives hints on a possible elongated shape of the core of the disclination line 
in the presence of the magnetic field. The extension of the defect core {i.e. the actual value 
of Xo) can only be determined when including the core energy into the investigation. This 
will be performed in the following section. 

We now proceed by inserting the third order polynomial expansion (|T2p into the equation 
(^). By comparison of the coefficients for the first two powers in the y coordinate {i.e. y^, y^) 
we obtain two ordinary differential equations for the unknown expansion coefficients ^o{x) 
and C{x). It is convenient to change to a set of dimensionless variables by measuring all 
length scales in units of yo, 

x = yox, ^ = 2$o, T = ylC. (15) 

We also introduce the notation 

V 16 y 

The equations for \E' and F have the following form 

^ (1 + F cos ^) ^" + (1 - cos ^) -AT)-K Q*' 2 - (r - sin ^ 
+2K{r -'^') sin^ + 2Z(r-^)^' cos^ 

-K (r - ^) ^' sin $ - - sin ^ = 0, (16) 

and 

- (1 + cos ^) (r" - ^") - - (r - ^) ^" sin ^ + 3 (1 - Kcos ^) T 

2 2 

+K (r - ^) (^ - 4r) sin * - ^ Z (r - ^') - 2 (r - ^) (^ - 4r)) sin ^ 

+K {^' - AT') sin ^ - Z (r - ^) Q ^'^ - (r - ^)^) cos ^ 
+2K(r-^)(r'-^') cos^-2Z^'(r-^)2 sln^ 

-K^'(r-^)2 cos$ + 2Z (^(r - ^) (r - ^') + ^ ^' - 4r)^ cos^ 

-K (^{T - ^) (r' - + ^ ^' - 4r)) sin ^ - ^ (r - ^) COS ^ = 0. (17) 

In (|T6|) and (|l^) ' denotes derivatives with respect to the dimensionless variable x. 

The set of ordinary differential equations ([l6|) and ([T7|) becomes much simpler for the one- 
constant approximation {K = 0). In this particular case the equations above are equivalent 
to the following ones 

^" = 8r-2^ + 2 sin^, (18) 
r" = 2r - 2^ + 2(r - ^) cos^ + 2 sin^. (19) 



Nevertheless, we shall analyse the set ([Tq) , ([171) . It turns out that a numerical solution and, 
if K is not too large, also an approximate analytical solution can be obtained. 
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According to ([I3D and (0) the boundary conditions (for positive solitons) now read (with 

^(x = 0) = 7r, ^(x = xo)=0, (20) 
r(x = 0) = -, r(x = Xo) = 0. (21) 



Eqs. (|TB|), (P7|), ([2DI), ( piD define a standard two-point boundary value problem. It 
can be solved numerically, for instance by a shooting method p!8|JT9| , pT|| . Satisfying the 
boundary conditions at x = 0, the ordinary differential equations (|T^), ( P7| ) are integrated 
numerically up to Xq. The integration constants are adapted iteratively in order to minimize 
the discrepancy between the numerical solution and the boundary conditions at Xq- For 
obtaining solutions for \E'(a;) and r(a;) we used a computer code from Numerical Recipes 



111 . Our calculations were performed for parameters corresponding to the liquid crystalline 
materials A^-(p-methoxybenzylidene)-p- buthylaniline (MBBA) and p-azoxyanisole (PAA) 
{see section IV). In these cases the numerical solution almost coincides with the approximate 
analytical solution presented below. 

An approximate analytical solution of (|16]), (|1^) can be achieved, which turns out to be 
quite accurate, as being revealed by a comparison with the numerical solutions. We start 
from the observation that the free energy density of the defect core, which corresponds to a 
disordered phase, is much higher than the typical elastic energy of the nematic. Therefore 
we expect that the core size xq is small in comparison with the half-width of the Neel wall 
?/o, i.e. Xq = Xq/uq <^ 1. Furthermore, ^ and F change by a finite amount on the interval 
[0,Xo], namely by vr or |, respectively. Therefore, the derivatives \Ef', F' are of the order 
Tr/xQ. They are much larger than ^> and F. One cannot a priori exclude that also the second 
order derivatives are large, of the order tt/xq. The approximation consists of keeping in (p!6[), 



(0) the leading terms only. Then ([l6|) reduces to 

(l+;^cos^)^"-^Z^'2 sin^ = 0, (22) 

which can immediately be integrated yielding 

(1 cos ^)^/2^' = const. (23) 

In the same approximation ( |T7|) simplifies to 

{l + K cos (F" - m") -Km' (F' - m') sin ^ 

(F - m) COS ^ - Z (F - m) m" sin ^ = 0. (24) 

In addition to the smallness of Xq ons can also exploit the fact that the elastic anisotropy 
K can be rather small. For example, for MBBA its value is 0.11, while for PAA it is 
0.42. Moreover, in (|2^) , (|^) K is multiplied by the sinus or cosinus of - this effectively 
diminishes the significance of the terms proportional to K even further. Therefore, it is 
natural to look for solutions of (|2^) , (^^ in the form of an expansion into powers of K. Up 
to first order in K we obtain 



^ = TT ( 1 - — 1 --K sin 



Xo/ 2 



vr (1-^ 

Xo 



+ 0{k'), (25) 
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and 



r = - (1 - ^ ) - 2 if sin 

4 V Xo 



X 

TT ( 1 - — 



+ ^^K (1-^1 (1 + cos 



( 1 ^ 

vr I 1 — — 

Xo 



) + OiK^). (26) 



A comparison with the numerical solutions of equations (p2D , (p^ shows that the functions 
in (p5|), (p6D yield very good approximations up to K = 0.6. 

IV. LENGTH OF THE DEFECT CORE 

Up to this stage, the length of the planar soliton xq (or, equivalently, xq = xo/yo) is 
unknown. We fix it by minimizing the total free energy, which includes the elastic and 
magnetic energy of the nematic phase as well as the energy of the defect core. 

Let us first calculate the total nematic energy -Fnem (per unit length along the z axis) 
for the soliton extending over the rectangle < x < xo, —yo < 2/ < l/o? which contains 
the core of the defect at ?/ = 0. Outside this rectangle there is the planar Neel wall at 
X < 0, —yo < y < yo, and the homogeneous director orientation parallel to the external 
magnetic field along the three remaining sides of the rectangle. Therefore the rectangle 
contains the total elastic and magnetic energy of the distorted nematic due to the presence 
of the defect. It is given by the integral 

F,en. = 2 / dx / dyj^^,^[<!>{x,y)], (27) 
Jo Jo 

where J-'nem is given by (^. For the tilt angle $(x, y) we use the approximate solution 
according to (|12|)> (HH)) (p5|) and (p6|). The integrals in (p7|) can be calculated by help of a 
computer algebra system {e.g., Maple). The result has the following form 



Fnera = - (Kn + K^z) 
2 



^ - 0.72 Xo + K f ^ + 0.93 - 11.69xo 

Xo V Xo 



(28) 



The terms proportional to 1/xo stem from elastic energy terms in J^-a_era proportional to 
Due to these terms the elastic energy of a defect with a point-like core would be infinite, 
because for such a defect xo = 0. 

The expression (p8|) would suggest that xq should be as large as possible - then Fnem 
would be minimal. In fact, this is not the case, due to the very large free energy stored in 
the defect core, where local transitions into disordered phases may occur. Let us perform an 
estimate of this core energy. (Again, it is understood that we consider the energy per unit 
length along the z axis.) The energy of the core is due to large gradients in the orientational 
order, which appear on a molecular length scale. Therefore it cannot be expressed in terms 
of the mesoscopic director field, but it is related to the molecular interaction potential across 
the discontinuity in the tilt angle on the segment < x < Xq of the x axis. This molecular 
interaction energy is small at the beginning and at the end of the core where the molecules 
on both sides of the segment are almost parallel. However, inside the core the molecules can 
even be perpendicular to each other. In this latter case the discontinuity of the tilt angle is 



8 



equal to |, and the separation of centers of mass of the molecules is of the order ^ where 
o"o denotes the molecular length. In the present paper we shall be satisfied with a rough 
estimate obtained by assuming that the core energy density is given by (0), when all terms 
are neg lected except for \ {Ku + K^^) with ^ f/;^. The width of the core is taken 
to be of the order This yields an estimate for the total energy of the core 



^2 



Feo.e ~ - (Kn + Ks^) ^1x0^ = -^ {Kn + i^ss) - xq. (29) 
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The total energy of the planar soliton is then F = F^em + -^corc- We now insert pq), (pOD 
and then minimize F with respect to the reduced core length xq. It is easy to find out that 
the Xo corresponding to the minimum total energy is given by 

° 1.7Ayo/ao-0.72-ll.Q9K' ^ ' 

Let us compute xq for particular nematic materials. For A^-(p-methoxybenzylidene)- 
p-buthylaniline (MBBA) at 25°C |3 the elastic constants are Ku = 6.0 ■ lO'^^ ^ and 



= 7.5 ■ 10~^^ N. The magnetic anisotropy is Ho^X = 9-7 ' 10^*^ Vs/Am, the molecular 
length o"o = 30 A. The magnetic field strength Hq is chosen 500 Oersted, according to a 
magnetic flux density Bq = ^iqHq = 0.05 T. Then, the elastic anisotropy is = 0.11. 
Equation yields yo = 3900 A. Finally, ^ 0.0027, and xq ^ 202 A. 

For p-azoxyanisole (PAA) at 120°C ||T^ the elastic constants are Ku = 7.0 ■ 10^^^ N and 



= 17.0 ■ 10 N. The magnetic anisotropy is /io Ax = 12.1 ■ 10 ^ Vs/Am, the molecular 
length (To = 20 A. The magnetic field strength Hq is again chosen 500 Oersted. The elastic 
anisotropy is = 0.42, yo = 6850 A, and finally xg ^ 0.0011, xq ^ 229 A. 

We notice that in both examples Xq is rather small, indeed. This is consistent with the 
assumption leading to the approximate solutions (pSf), p6|). The resulting physical length 



of the core xq is relatively large and it probably could be seen in appropriate experiments. 
The dependence of the nematic, core and total energies on the reduced length of the defect 
core is plotted in Figs. 3 (MBBA) and 4 (PAA). 

With the determination of the reduced core length xq the calculation of the director field 
for the planar soliton is completed. The tilt angle field is shown in Figs. 5 (positive soliton) 
and 6 (negative soliton). The core line at ?/ = is clearly visible by the jump of the tilt 
angle. However, this picture is somewhat misleading, because it does not take into account 
that the director is an object without arrowhead. For instance, at x = 0, ^ = 0± there is a 
jump by TT which in fact means an orient at ional change of zero angle, exactly the same as for 
X = Xo = 0.052, y = 0±. As already stated in the previous section, due to the periodocity of 
TT for tilt angle changes the largest orientational jump occurs for Xd = 0.025, y^ = 0±, where 
the tilt angle is ±|. This point can be defined as the center of the core which is related to 
the original disclination line (in three dimensions). 

These particular features become obvious from a lattice visualization of the director field, 
which is presented in Figs. 7 (positive soliton) and 8 (negative soliton). (In these figures the 
X and y dimensions are not proportionally scaled.) Rods of unitary length placed on the 
sites of a rectangular lattice indicate the local orientation. The dashed line means the defect 
core and the small circle marks the center of the core, according to the previous discussion. 
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V. REMARKS 



The positive and negative planar soliton are distinguished by the boundary conditions 
( pOD and (pID, but not by the field equations. Therefore they cover exactly the same 
area, although their energy content is slightly different. The tilt angle field for the 
negative soliton is obtained from the positive soliton solution y) pre- 

sented above by a sign inversion of the expansion coefficients: $o (x) = —^o{x) and 
C-{x) = -C{x). 



2. Expression (pOj) reveals a rather interesting dependence of the reduced core length 
Xq on the magnetic field Hq which enters through yo {see ([TT])). For weak magnetic 
fields we have large yo, hence xq is small and it tends to zero when the magnetic field 
vanishes. However, the physical length of the core is equal to xq = yox and it increases 
as ^^7= when the magnetic field decreases. The physical reason is that for a weaker 
magnetic field the distance over which the director field can be reorientated by a given 
angle is larger. In the case of planar solitons the required reorientation is such that 
the tilt angle changes from $Neei(l/) towards zero. Of course in this limit the width 
of the Neel domain wall (|lT]) also increases, as -j^, hence faster. On the other hand, 
with increasing magnetic field yo decreases and Xq increases. From inserting ([TI|) into 
( pO] ) it is noticed that formally there is a finite critical value for the magnetic field at 
which Xq becomes infinite. For MBBA the dependence of the reduced and physical 
core length xq and Xq on the reduced magnetic field h = if [Oersted] /500 is 



To = 0,55,/ " ,, f^.^M^ (31) 
"112,8 -ft- 202 ./ft (112.8 -A) ^ ' 



which yields a critical reduced magnetic field = 112.8, corresponding to a critical 
fiux density of 5.64 T. (Analogous calculations for PAA give a somewhat smaller 
critical fiux density of 5.3 T.) However, it should be remembered that Eq. (|30|) has 
been derived under the assumption that is small, so it may become wrong well 
before reaching the critical value of the magnetic field. 

To conclude, we found an approximate solution for the director configuration in planar 
solitons in nematics. It is continuous everywhere apart from a strip of finite width. This 
points out the possibility of an elongated shape of the defect core in disclination lines due 
to an external magnetic field. 
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FIGURES 



FIG. 1. Geometry and coordinates for a positive planar soliton in a nematic liquid crystal. 

FIG. 2. Geometry and coordinates for a negative planar soliton in a nematic liquid crystal. 

FIG. 3. Free energy (per unit length) F' vs. core length xq for MBBA at 25°C. Both 
quantities are in reduced units (dimensionless). F' = F/Kscv, xq = xo/yo, with units 
i^av = \ {Ku + K-is) = 6.75 • 10-12 N, yo = 0.39 //m. Dashed line: analytical solution for the 
energy of the nematic phase; rhombs: numerical solution for the energy of the nematic phase; 
dotted line: analytical solution for the energy of the defect core; crosses: numerical solution for 
the energy of the defect core; solid line: total energy. 

FIG. 4. Free energy (per unit length) F' vs. core length xq for PAA at 120° C. Both 
quantities are in reduced units (dimensionless). F' = F/Kgy, xo = xo/yo, with units 
K^v = \ {Kn + Kss) = 12 • lO-^^ N, yo = 0.68 jj^m. Dashed line: analytical solution for the 
energy of the nematic phase; rhombs: numerical solution for the energy of the nematic phase; 
dotted line: analytical solution for the energy of the defect core; crosses: numerical solution for 
the energy of the defect core; solid line: total energy. 

FIG. 5. Tilt angle field for a positive planar soliton in MBBA at 25°. Spatial coordinates in 
reduced units, x = x/yo, y = y/yo (yo = 0-39 fxm). 

FIG. 6. Tilt angle field for a negative planar soliton in MBBA at 25°. Spatial coordinates in 
reduced units, x = x/yo, y = y/yo {yo = 0.39 fim). 

FIG. 7. Lattice visualization of the director configuration for a positive planar soliton. 

FIG. 8. Lattice visualization of the director configuration for a negative planar soliton. 
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